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Abstract 

Energy saving is becoming an important issue in the design and use 
of computer networks. In this work we propose a problem that considers 
the use of rate adaptation as the energy saving strategy in networks. The 
problem is modeled as an integral demand-routing problem in a network 
with discrete cost functions at the links. The discreteness of the cost 
function comes from the different states (bandwidths) at which links can 
operate and, in particular, from the energy consumed at each state. This 
in its turn leads to the non-convexity of the cost function, and thus adds 
complexity to solve this problem. We formulate this routing problem as 
an integer program, and we show that the general case of this problem is 
NP-hard, and even hard to approximate. For the special case when the 
step ratio of the cost function is bounded, we show that effective approx- 
imations can be obtained. Our main algorithm executes two processes in 
sequence: relaxation and rounding. The relaxation process eliminates the 
non-convexity of the cost function, so that the problem is transformed 
into a fractional convex program solvable in polynomial time. After that, 
a randomized rounding process is used to get a feasible solution for the 
original problem. This algorithm provides a constant approximation ratio 
for uniform demands and an approximation ratio of 0(log''~^ d) for non- 
uniform demands, where /3 is a constant and d is the largest demand. 
Keywords: energy saving, network optimization, network routing, ap- 
proximation 



1 Introduction 

Energy-aware computing has recently become a hot research topic. The increas- 
ingly widespread use of Internet and the sprouting of data centers are having 
a dramatic impact on the global energy consumption. The energy consumed 
comes from the aggregate power used by many devices (CPUs, hubs, switches, 
routers). As shown in [?], up to 45% of the total energy is consumed by com- 
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puting resources such as CPUs and storage systems. There has been significant 
amount of work done on energy management schemes for these largest energy 
consumers, which was started by the work on power consumption minimiza- 
tion in processors done by Yao et al. [?]. In parallel with energy saving in 
computing devices, the energy consumption of network devices has also been a 
fundamental concern in wired networks like Internet or data center networks. 
Gupta and Singh [?] were among the first to identify energy consumption of 
general networks as an important issue. Kurp et al. [?] then showed that there 
is significant room for saving energy in current networks in general. The main 
reason is that networks arc always designed with a significant level of rcdim- 
dancy and over-provisioning, in order to guarantee QoS, and to tolerate peak 
load and traffic variations. However, since networks usually carry traffic that is 
only a small fraction of the peak, a significant portion of the energy consumed 
is wasted. This is particularly representative in data center networks, where 
the connectivity redundancy is very heavy and the traffic load is dramatically 
fluctuant with certain patterns over time. Ideally, the energy consumed in a 
network should be proportional to the traffic load it carries. 

In this paper we consider a centralized energy saving problem for networks. 
In principle, this problem has applications in most networks, and can partic- 
ularly be used in data center networks or other software-defined networks. In 
the model assumed, we use an energy saving strategy called rate adaptation , 
which assumes that the network links have a discrete set of operational states 
(bandwidths) , each with its corresponding power consumption. By using the 
rate adaptation strategy, network devices choose an appropriate speed according 
to their current traffic load. This model is different from the usual speed scal- 
ing model (see below) that assumes a continuous power function, and is more 
realistic. However, the optimization problem we need to solve in a model with 
rate adaptation is more complex compared with using these other models. Gu- 
naratne et al. [?] proposed a network model in which links use rate adaptation, 
that they called adaptive link rate (ALR). We will present a formal description 
of our network model, that uses rate adaptation as power saving strategy, and 
of the problem of providing route assignment to demands from a global view of 
the network. Actually, this model turns out to be a network routing problem 
with a discrete cost function. We show by analysis that in the general case, this 
problem is hard to approximate. Given some restriction on the cost function, 
we provide an efficient approximation algorithm for it. 

1.1 Related Work 

Work on energy efficiency have mostly focused on two strategies to save energy: 
speed scaling and powering down. Under speed scaling, it is assumed that the 
power consumed by a device working at speed s has the form of P = s'^ , where 
^ is a constant between one and three. This comes from the well known cube- 
root rule, which states that the speed is approximately the cube root of the 
power consumed. Speed scaling has been widely used for optimizing the energy 
consumption of computing resources, especially the processor (see [?, ?, ?, ?]). 
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This strategy has also been used for network devices. However, most studies 
in this category focused on a single device. At the same time, speed scaling 
as described does not model realistically network devices and is hard to be 
applied in practice. Another approach to save energy is achieved by powering 
down the devices while they are idle. Andrews et al. [?] considered that network 
elements operate only in the full-rate active mode or the zero-rate sleeping mode. 
Nedevschi et al. explored both speed scaling and power down to reduce global 
network energy consumption. ElasticTree [?] is a centralized energy controlling 
system for data center networks which turns off some of the routers or switches 
and then yields energy-efficient routes. 

The general network routing problem is described as follows. We are given 
a set of traffic demands and want to unsplitably route them over a transmission 
network. The total traffic on link e incurs a cost which is defined by a cost 
function fe{xe)- Our objective is to find routes for all demands such that the 
total incurred cost fe{xe) is minimized. There has been significant work 
on this general network routing problem which is also called the minimum-cost 
multi-commodity flow (MCMCF) problem. Note that the complexity of this 
problem depends on the cost function defined on each edge. In the simplest 
case where the cost function fe{xe) is linear with the carried traffic .Tg, this 
problem turns out to be the shortest path problem which is polynomial-time 
solvable. Moreover, if we choose fe{xe) as a subadditive continuous function 
which has the property of economies of scale, the problem becomes the Buy- 
at-Bulk (BAB) problem which has been well studied. Awerbuch and Azar [?] 
provided an O(log^n) randomized approximation algorithm for this problem. 
Andrews [?] showed that for any constant 7 > 0, there is no 0{\og^~'' n)- 
approximation algorithm for uniform BAB (all edges have a uniform cost func- 
tion), and there is no O{\og^ ' n)-approximation algorithm for the non-uniform 
version, unless NP e ZPTIME{nP°^y'"s Here n is the size of the network. 
However, for general cases with superadditive continuous cost functions, achiev- 
ing any finite ratio approximation is NP-hard [?]. Andrews et al. [?] also studied 
an energy-aware scheduling and routing problem for wired networks under the 
speed scaling model, presenting a constant approximation algorithm for uniform 
demands. In [?, ?] the authors studied routing under a special case of speed 
scaling with a startup cost, and presented a polylogarithmic approximation al- 
gorithm. The most significant difference between the work mentioned above and 
our work is that here we focus on the network routing problem with discrete 
cost functions rather than continuous ones. 

1.2 Our Results 

Based on the general network routing problem, our problem can be described as 
follows. We want to route a set of traffic demands over a transmission network. 
Each link (which abstracts both the transmission link and the corresponding 
ports of the end devices connected by this link) can operate at one of a given 
set of speed states Ri < ■ ■ ■ < Rm- For each speed state Ri, a fixed amount of 
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cost Ei is defined to represent the power consumption at this speed. To carry 
ah traffic flows, each hnk e chooses a speed state Ri at which to operate so that 
Xe S {Ri-i,Ri]. If Xe < Ri, then i?i is chosen. We also assume that Xe < Rm 
for any routing considered. Consequently, each link will incur a cost and the 
total cost of the network can be defined as Ef, where Ef is the cost of link e 
while choosing a speed Ri. Equally we can use a cost function fe{xe) to describe 
the cost on each edge. This cost function can be obtained by setting the value 
of fe{xe) to Ei if Xe & {Ri-i,Ri]. Thus the total cost can be represented as 



Xle/e(^e)- It is trivial that function fe{xe) has a discrete formulation and is 
indeed a step function. For simplicity, we focus on the uniform case, where all 
the links consume energy in the same way and thus have the same cost function 



We aim to solve the minimum-energy routing with this discrete cost function 
in this paper. In Section [21 we give the formalized description of the model and 
prove that in general we cannot get any finite ratio polynomial-time approxi- 
mation algorithm for it if P^^NP. In Section |31 we consider a special case where 
/(I) and the ratio between any two adjacent steps of the cost function /(■) is 
bounded by a constant. Given this assumption, we show that efficient approxi- 
mation algorithms can be achieved. The approximation we provide in this paper 
is established on a two-step relaxation and rounding process. The first relax- 
ation is made on the cost function to eliminate the discreteness, which is done 
by a special interpolation method which transforms the discrete function /(•) 
into a continuous one g{-) while introducing a bounded error. Then, by relaxing 
the integral indicator variable, the induced problem turns out to be a convex 
program and can be solved in polynomial time. Hereafter, a two-step rounding 
technique is used to round infeasible solutions to the feasible region. We show 
by analysis in SectionUthat our approach can achieve a constant approximation 
for the cases with uniform demands and an 0(log^~^ d)-approximation to the 
cases with non-uniform demands, where d is the maximum one among all traffic 
demands. Finally in Section [SI we draw conclusions and provide some future 
work and open problems derived from this work. 



We study the energy-efficient routing problem under a general case where the 
optimizer takes a traffic matrix and a network topology as inputs and returns 
the network configuration. We first consider a simple case with unit demands 
where all the demands have a flow amount of one integer unit, and then extend 
the results to the uniform and non-uniform cases, where we consider that all 
the traffic demands have the same flow amount of some integer units and have 
different flow amounts, respectively. 

The problem is described as follows. Assume we are given an undirected 
graph G = {V,E) and a set of traffic demands D = {di, ^2, rffc}, where the 
i*'' demand, 1 < i < fc, requests di units of bandwidth provisioned between a 
source node Si and a sink node U. Unless otherwise stated, in the following we 
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2 The Model 
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assume unit demands, i.e., di ^ 1. We also assume that links represent abstract 
resources (links with their switch ports), and each link can operate at one of 
a constant number of different rates Ri < R2 < ... < Rm- Here we assume 
that Ri > 1 because a link will carry a flow of at least one unit. For energy- 
saving consideration as we mentioned in the previous section, it is reasonable to 
have many different rates in future network devices. Each rate Ri, 1 < i < m, 
has a cost of fe{Ri) which is the power consumption. Our goal is to route 
every demand in an unsplittable fashion with the objective of minimizing the 
total cost X^g/e(-Re)) where Re e {Ri\i S [l,...,m]}. Note that unsplittable 
routing is important in networks in order to avoid packet reordering. Wc only 
consider the case in which all links have the same cost function /(•). We call 
this energy-efficient routing problem with discrete cost functions as rate adaptive 
energy- efficient routing problem. 

Not surprisingly, the rate adaptive energy-efficient routing problem is NP- 
hard, due to the fact that the cost function is discrete. Furthermore, we show 
that, in general, it cannot even be approximated. This is shown by the following 
theorem. 

Theorem 2.1. For any constant p > 1, there is no polynomial-time p- approximation 
algorithm for the rate adaptive energy-efficient routing problem, unless P=NP. 

Proof. We show here that a polynomial-time p-approximation algorithm ^ 
could be used to solve the edge-disjoint paths (EDP) problem, which is NP- 
hard, proving the claim. Let us assume the algorithm exists. Let us consider 
an instance of the EDP problem, given as a network with w links and a collec- 
tion of k pairs source-sink, {sj,tj). We transform it into an instance of the rate 
adaptive energy-efficient routing problem as follows. In this instance, the net- 
work, the number of demands, and the source Sj and sink tj of each demand j 
are the same as in the EDP instance. Also, there are only two rates Ri and R2, 
R2 > k ■ Ri, and we choose any cost function such that /(i?2) > p ■ w ■ f{Ri). 
Finally, we set each demand value dj = Ri. This transformation is trivially 
polynomial in time. 

We use the algorithm ^ to find a solution that p-approximates the optimal 
solution of the rate adaptive energy-efficient routing problem obtained by the 
transformation. (Observe that since R2 > k ■ Ri, a solution always exist.) Let 
C be the cost of the solution found. We claim that C < p-w- f{Ri) if and only 
if there are edge disjoint paths for all the demands. If this claim holds, then ^ 
can be used to solve EDP in polynomial time and hence P = NP. 

The claim can be proven as follows. (=J>) Assume that C < p ■ w ■ f{Ri). 
Then, the paths found by ^ are edge disjoint. Otherwise, at least one of the 
links would be used by at least two demands, the rate of that link would be 
R2, and C > /(-R2) > p ■ w ■ f{Ri). {<=) Assume edge disjoint paths exist for 
all demands. Then, there is a way to route all the demands so that each link 
is used by no more than one demand, and can be set to rate The cost of 
every edge in this solution is at most f{Ri), and the cost of the solution is hence 
C < w ■ f{Ri). Trivially, the optimal solution also has cost C* < w ■ f{Ri). 



5 



Since ^ is a p-approximation algorithm, the cost of the solution it outputs must 
satisfy C<p-C* <p-w/(i?i). □ 



Actually Lemma 2 in [?] is a special case of this theorem, in which in the proof 
-Ri = 1 and f{Ri) = 0. The proof of the theorem suggests that we cannot get 
any finite-ratio polynomial-time approximation algorithm for the rate adaptive 
energy-efficient routing problem unless we bound the ratio of costs of different 
rates and /(l)/^ (we call it as step ratio). In the rest of this 

work, we will considcir that this ratio is bounded by a constant cr. Note that this 
is quite sensible because the energy consumption of devices is always bounded. 
But, even in this special case this problem is NP-hard. The good news are 
that the problem with constant step ratio can be approximated by a constant 
approximation ratio. We will show the details in the following sections. 

Formally, we formulate the described routing problem into an integer pro- 
gram, that can be seen as (Pi) below. The binary variable jji^e indicates whether 
demand i uses link e, while Xe is the total load on e. Flow conservation means 
that for each demand i the source Si generates one unit of flow, the sink absorbs 
one unit of flow, and for the other vertices the incoming and outgoing flows of 
demand i are the same. Observe that for Xe < Zg, f{xe) = f{ze)- This results 
in the discrete property of the cost function /(•). More precisely, f{x) is a non- 
decreasing step function of x, where x is the speed of each link. In practice, 
cost functions for network resources might different, but here we focus on the 
case with only a uniform cost function. There is no doubt that solving (Pi) 
is NP-hard for the — 1 constraint on variable yi_e- Since solving our network 
routing problem is NP-hard, we have no hope on finding the optimal solution. 

(Pi) min ^f{ze) 

e 

subject to 

i 

^ Zg 

Ze G {-Rl, Rm} 

yi,e e {0, 1} 

Ui^e '■ flow conservation 

3 Approximation 

In this section, we present our main algorithm, which is devised to approximate 
the optimal solution of (Pi). As indicated above, the complexity of the rate 
adaptive energy- efficient routing problem comes from the non-convexity of the 
cost function /(■) and the integral constraint on j/i^g. Intuitively, if we remove 
these constraints and transform function /(•) into a convex one, the induced 
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problem could be solved. From this intuition, we propose an approximation 
algorithm for this problem. By a two-step relaxation and rounding procedure, 
the algorithm achieves a constant approximation. 

In order to get an approximated solution of (Pi), a two-step relaxation and 
rounding procedure is introduced. The first relaxation is made on the cost func- 
tion /(■), to remove the discreteness of /(•). We use a particular interpolation 
method to transform the discrete cost function of the original program into a 
continuous convex one. It makes the program to be simpler while introducing a 
bounded error. Additionally, we relax the binary variable yt^e to be real. Then, 
the induced problem is a convex program which has been proved to be optimally 
solvable in polynomial time. Now we describe our algorithm in detail. 

3.1 RelcLxing the Cost Function 

The first relaxation is made on the cost function. We use a special interpolation 
method to simplify (Pi) by replacing the step function /(•) with a continuous 
function .g(-). 

Before applying the interpolation, it is important to decide the form of the 
function g(-) which we aim to obtain. Due to the well known cube-root rule, 
it has been suggested that most network devices consume energy in a super- 
additive manner ([?, ?]). That is, doubling the speed more than doubles the 
energy consumption. In particular the energy curve is often modeled by a poly- 
nomial function g(x) = fix^ , where fi and /3 arc constants associated to the 
network elements. More precisely, the parameter /3 in the ordinary form of en- 
ergy consumption has been usually assumed to be in the interval (1,3] [?]. The 
objective here is to transform a step cost function /(•) into a function in the 
form of g{x) = jix^ . Although, as mentioned, /3 will be typically larger than 
1, and hence g{-) will be a convex function, the proposed interpolation method 
does not impose any restriction on this. 

Now we discuss how to apply the transformation from a step function to a 
continuous (convex) one. A common approach is to use the midpoints of all the 
steps as discrete values and fit by mean squares. This approach will not give 
good results if the steps of the function have different lengths, because there 
might be many steps with small lengths dominating the interpolation. Another 
popular method is to perform an interpolation on a set of points which is ob- 
tained by sampling the original function. Unfortunately, using this technique 
the error of the interpolation depends on the sampling method we choose, and 
is always hard to estimate. In order to contain the interpolation error, we devise 
a new interpolation method which based on minimizing the difFcrcnce between 
the two function. Without depending on some other parameters, the method 
works well especially for the fitting of step functions. 

Consider the original function f{x), and the one to be fitted g{x), as previ- 
ously mentioned. f{x) is defined as follows. 
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2/1, l<2:<i?i, 

(1) 

2/rn , Rm~l <^ 3^ ^ Rm^ 

where in our case yi = f{Ri), {1 < i < m) is the cost of each state and Ri, Ri+i 
(1 < z < to) represent the lower and upper boundaries of the speed for each 
state. We aim to fit g{x) to f{x). 

The formula that has to be minimized can be represented as 

xell,R„^] [ g{x) f{x) J 

We only consider x G because /(O) = ^(O) and in any feasible integral 

solution, X ^ (0, 1). Since g{x) is not a linear function, this minimization 
problem is hard to solve. We then consider an alternative. We use the fact 
that the formula becomes linear with parameters fi and /3 when a logarithmic 
transformation is applied. The fitting function g(x) is clearly linear under a 
logarithmic transformation. Note that 

log(5(a;)) = log M + /? log X. (3) 
Then it is equivalent if we minimize the following formula 

max I log /(a;)- log I 

= max I log/(a;) — (log^ + /31oga;)|. 
Ee[i,i?m] 

However, the above formula of G' is still unable to tackle, because the absolute 
operation is hard to handle. Here we propose to use an alternative which is to 
minimize the integral of the square of the difference between the two functions. 
Let us define Vi = logy^, Wi = log Ri for i G [1,to] and wq = logl = 0, and 
/z' = log fi. Then the alternative formula that we will in fact use is 

H{ti', f3)^Y. - + M?dw. (5) 

i=i "'"'■-I 

And now Eq. ([5]) is to be minimized with respect to the parameters of the general 
quadratic equation. We choose fj,' and /? such that the first partial derivatives 
of H{iJ.',l3) are equal to zero and its second derivatives are positive. 



dH 



'2[vi ~ fJ.' — /3w]dw = 0, 
-2w[vi — ii' ~ /3w]dw = 0. 



(6) 



8 



It is obvious that the second derivatives are all positive. By solving Eq. we 
get the values of parameters fi' and /3, and from /i' we obtain /j,. Consequently, 
the objective function g{x) of the interpolation can be determined. Once the 
function g{x) is obtained, the optimization problem can be rewritten as follows. 

(P2) min ^g{xe) 

e 

subject to 

Xe=^ Vt.e Ve 

i 

yi,e e {0, 1} Vi, e 

Ui^e '■ flow conservation 

The problem now turns to be an integer program with a convex objective func- 
tion. Solving (P2) is stiU NP-hard. 

3.2 Relaxing the Binary Constraint 

In this section, we show that by relaxing the binary constraint on y^ e, P2 can be 
efficiently solved. Although this brings some accuracy loss, it makes the problem 
solvable. It is feasible to do this kind of relaxation in our case, with some small 
modifications. One key observation is the convexity of the objective function in 
P2. By combining the linear constraints and relaxing the binary variable yi^e 
from {0, 1} to [0, 1], P2 turns to be a regular convex program and can be solved 
in polynomial time by using algorithms for convex programming. We denote 
the fractional solution obtained by convex programming as y* ^. However the 
problem is that in the fractional solution y* ^, the trafflc flow of a demand can 
be splitted over multiple paths, which is not feasible in our case. We address 
this below. 

3.3 Two-step Rounding 

In this section, we introduce a two-step rounding technique to transform the 
fractional solution into a feasible one. First we round the fractional routing 
solution to an integral one, and then we determine the link rates from the 
routes. 

For the routing path rounding, the key point is how to transform the frac- 
tional routes into a solution in which there is only one path for each traffic 
flow. In the fractional solution, the flow i carried by link e is represented as y* ^, 
where y*^ G [0, 1]. Our goal is to get a path pi for each demand i which satisfies 
{e|e S Pi} C {e\y* ^ > 0}, to route flow i. We use the Raghavan-Thompson 
randomized rounding method to complete this transformation. 

The overall rounding process is described below. Once the optimal fractional 
solution y* ^ has been found, the flows assigned to the links is mapped to paths 
as follows. For each demand i, first we generate a sub-graph Gi defined by links 
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e where y* ^ > 0. Then, we extract a simple path p connecting the source and 
destination nodes. The hnk e € p with the smaUest weight ^ is called the 
bottleneck link. This y* ^ is selected as the weight of this path, denoted as Wp. 
Hereafter the weight y*\^ of each link in path p is decreased by Wp. The above 
procedure will be repeated until that for all e we have y*^ = 0. By the flow 
conservation constraint, we can state that 

Proposition 3.1. By repeating the path extraction process on Gi, a state in 
which ye G Gi, we have y*e = reached. 

Consequently, we define this state as the termination of the path extraction 
process. As a result we get a collection of paths {p} and for each of them we 
have a weight Wp. We randomly select one path from {p} for current demand 
i using the path weights as the selection probabilities. After this rounding, 
there will be only one single path for each demand. We denote this solution as 
vie e {0, 1}. 

After the routes for the demands have been chosen, the state of every link 
must be determined. We select the speed of each link via the following rounding 
procedure. First, we compute the carried traffic Xe = J^iVle, where yi^e is the 
amount of demand i that traverses link e after the rounding. Then for each link, 
we search the collection of possible operational speeds and choose the minimal 
Se which can support the carried traffic. Formally we have 

Se = mm{Ri\{i G [l,m]) A {x'e < Ri)}. (7) 

Determining the link states and routing all the traffic demands result in an 
approximated solution for the rate adaptive energy-efficient routing problem 
with discrete cost functions. 



4 Algorithm Analysis 

We evaluate our method in this section. As it is designed to approximate the 
optimal, it is essential to derive an approximation ratio as the main accuracy 
criterion. Some accuracy loss has been introduced in both relaxation and round- 
ing processes, so we study them one by one. First, we show a bound on the 
error introduced by the interpolation. Second, we bound the error by rounding. 
Then, by combining these results together, we draw our main conclusions. 

4.1 Bound on the Interpolation Error 

Interpolation is proposed to approximate the objective function, so it is impor- 
tant to bound the error introduced. During the interpolation process, the error 
comes from the gap between the original function f{x) and the fitted function 
g{x). We define this gap as follows: 

gap= max S — tT' TTT r • (8) 
xe[i,Rm] { g{x) f{x) j 
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Here we also consider x G as we did in previous sections. We will use 

this gap definition as the interpolation error, reason for which we present the 
following theorem. Here we assume that g{x) intersects with f{x) in each step 
of f{x), which is reasonable because we assumed that g{x) could well describe 
/(x) under the cube root rule. If not, there could be a big difference in trend 
between f{x) and g{x), and we cannot obtain any bounds. In other words, this 
would mean that the cube root rule does not apply to this network. 

Theorem 4.1. Given function f{x), if f{x) satisfies yi/yi-i < a where a > 1, 
then in interval [l,Rm], the interpolation error satisfies gap G [^pj-,iy9], where 
(f = max{cr,/(l)/^}. 

Proof. First we show the lower bound. Assuming f{x) and g(x) intersecting 
at each gap, and given the form of g{-), means that we can focus on values of 
X S {1, Ri, R2, Rm}- Hence we consider a particular point x from the set 
{1, i?2, -Rm} (see Fig. [T]). Assume that we can bound the gap and that 
we have a bound S. Let 7 — maxi{yi/yi_i}. 

Case 1: The fitted function g{x) is above the midpoint of the two values of 



the step function f{x) at x, as can be seen in Fig. 1(a) Therefore, at x, the gap 
is defined as gap — g{x)/ f{x). Suppose we are given S < ^^^X.. Then for 7 > 0, 
we have 

-^<^- (^) 

Then we can obtain 

7= — > (2(5-1). (10) 

Vi-i 

Notice that yi equals f{x) and here g{x) is not smaller than + . We 

have 



gix) > yi 



Vi - Vi-i 



2 (11) 

Syi-i = Sf{x). 

This results in g{x) > Sf(x), which is a contradiction to the assumption that 
g{x)/f{x) < S. So we have S > in this case. 

Case 2: The fitted function g{x) is below the midpoint of the two values of 



f{x) But x' ~ X + e where e is infinitesimal, as can be seen in Fig. 1(b) Here the 



gap is defined as gap — f{x')/g{x'). Suppose we are given 5 < We have 

2y,> ^iy^-l+y^). (12) 
Note that y, = f{x') and g{x') < y^^i + Then we can derive 

g{x ) < y.i-1 + 

Similarly to the first case, the inequality f{x') > Sg{x') results in a contradiction 
to the assumption that we have a bound. So we also get 5 > -^^^ in this case. 
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Ox Ox 

(a) g{x) > 2 (b) 9m < 2 

Figure 1: Two cases of error bounding proof 

Table 1: Notation for Bounding the Approximation Ratio 

Name Definition 

Cf optimal integral solution under /(■) 

Cg optimal integral solution under g(-) 

optimal fractional solution obtained by the program 
■f with a binary relaxation, which is fi^t) 

C* optimal fractional solution under g{-), which is gi^t) 

Cg solution obtained after the first rounding process, which is X]e3(^"e) 

C f solution obtained by our overall method 



As we have obtained the lower bound, now we discuss about the upper 
bound. It is straightforward that we cannot get a tighter upper bound rather 
than S for the error at any point x G [-Ri, i?2, ^m]- For the extreme case 
where the fitted function g{x) takes a value which is very close to yi at x, the 
gap can be represented by gap — g{x)/f{x) w yi/yi-i, as well as the error. 
However, when x G [1, this gap may be bounded by the case when a; = 1, 
where gap = Therefore, the interpolation error can be upper bounded 

hy (p = max{(7,/(l)//x}. 

Combining these results, we have that the interpolation error is in 'fi ■ 

This completes the proof. □ 
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4.2 Approximation Ratio Analysis 

In this section, we present our results for the approximation ratio. The notation 
that we are going to use are presented in Table [1] Given this notation, our 
objective is to bound the ratio between Cf and . We will assume that a and 
ip are bounded by constants and will show that the ratio is then bounded by a 
constant. 

Lemma 4.2 ([?]). For unit demands, randomized rounding on paths can obtain 
a X- approximation for the energy- efficient routing problem such that E[Cg] < 
XE[C*], where X is a constant. 

Since C* is obtained by the optimal fractional solution and it is always the 
lower bound of the optimal integral solution Cg, we directly have that E[Cg] < 
XE[C;] < XE[Cg]. 

Lemma 4.3. // the ratio between any two adjacent steps of cost function /(•) 
is bounded by a, then Cf < ipaCg. 

Proof. The result follows from Theorem 14. 11 in which it shows that the largest 
gap between g{X(.) and /(xg) is (f. Then, from the relation between Sg and 
(see Eq. ([7])), we have that /(se) < <^f[xe). By Theorem 14.11 /(sg) < (pag{xe) 
follows. Thus, the following result can be derived 

C/=^/(Sg)<(p(J^5(Xe)==:</JCr4- (14) 

e e 

This completes the proof. □ 

Theorem 4.4. For unit demands, the expected energy consumption Cf, is a 
p- approximation of E[Cf]. That is, E[Cf] < pE[Cf], where p is a constant that 
depends on a and ip. 

Proof. The expected energy consumption of the solution found is E[Cf] = 
E[J2e fi^e)]- From Lemma [4.31 we have that Cf < fcrCg and, hence, E[Cf] < 
<paE[Cg]. As it is shown in Lemma 14. 2[ there is a constant A such that 
E[Cg] < XE[C;]. 

To complete the proof, we observe from Theorem l4 . 1 1 that . for all x, g{x) /if < 
f{x). Then the optimal fractional solution under g{x) satisfies Cg/ip < Cj. 
Given that the optimal fractional solution under f{x) is always a lower bound 
of the optimal solution Cf, we also have that C* < (pCj < ifCf. The last 
inequality comes from the fact that the optimal fractional solution is always a 
lower bound for the optimal solution. Combining all these results together, we 
can obtain that 

E[Cf] < wE[Cg\ < f<TXE[C*g] < if'^aXE\Cf]. (15) 

Assume that p = Lp^aX. As a, ip and A are constants, we conclude that p 
will also be a constant. This completes the proof. □ 
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This result can be directly applied to uniform demands, where each traffic 
demand requests a bandwidth di = d. In this case the total flow on each edge 
will be d times that of the case with unit demands. This is the only difference 
between these two cases. 

Corollary 4.5. A p- approximation to the optimal integral solution in expecta- 
tion can be obtained for uniform demands. 

For non-uniform demands, Andrews et al. [?] has derived an approximation 
result on the randomized rounding on paths. We borrow this result and combine 
it with Theorem 14 . 1 1 and Lemma Then the following theorem holds. 

Theorem 4.6. For non-uniform demands, we can achieve a 0{log^^^d)-approximation 
for the rate adaptive energy- efficient routing problem, where d is the maximum 
traffic demand in D. 

5 Conclusion and Future Work 

In this paper, we investigate the energy-efficient network routing problem with 
discrete cost functions. Our contributions are mainly on the following results. 
For the rate adaptive energy-saving problem, we formulate it as an integer 
program and prove that the general case of this problem is NP-hard, and even 
hard to approximate. For a special but practical case of this problem, we provide 
an efficient approximation algorithm which is based on a two-step relaxation and 
rounding process. We show by analysis that the error incurred in each step of 
our algorithm can be bounded. By using our algorithm we obtain a constant 
approximation to the optimal for unit and uniform demands and an 0{log^^^d)- 
approximation for non-uniform demands. 

Observing that the network traffic comes from the end devices connected 
to the network, we have that it is important to combine the optimization of 
placing end devices with the network routing optimization. For instance in a 
datacenter environment, a good placing of virtual machines can significantly 
reduce the network congestion as well as the energy consumption of network 
devices. As a natural extension to our routing problem, this new problem can 
be described by an additional description to the original model saying that we 
can also determine the sources and destinations of all the demands as well as 
the routing paths. Actually this turns out to be a combination of a quadratic 
assignment problem and a network routing problem. In general, this problem 
is much harder than the one we discussed in this paper. To find effective ways 
to approximate it is still an open problem. 
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